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Finite thin disc models of four galaxies in the Ursa Major cluster: 
NGC3877, NGC3917, NGC3949 and NGC4010 
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ABSTRACT 

Finite thin disc models of four galaxies in the Ursa Major cluster are presented. The models 
are obtained by means of the Hunter method and the particular solutions are choosen in such 
a way that the circular velocities are adjusted very accurately to the observed rotation curves 
of some specific spiral galaxies. We present particular models for the four galaxies NGC3877, 
NGC3917, NGC3949 and NGC4010 with data taken from the recent paper by Verheijen & 
Sancici (2001). By integrating the corresponding surface mass densities, we obtain the total 
mass M. of these four galaxies, all of them being of the order of 1O 1O A^0. These obtained 
values for M. may be taken as a quite accurately estimative of the mass upper bound of these 
galaxies, since in the model was considered that all their mass was concentrated at the galactic 
disc. The models can be consider as a first approximation to the obtaining of quite realistic 
models of spiral galaxies. 
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1 INTRODUCTION 

Currently, the most accepted description of the composition of spi- 
ral galaxies is that a main part of its mass is concentrated in a thin 
disc, being the other constituents a spheroidal halo, a centra l bulge 
and, perhaps, a central black hole dBinnev & Tremainel2 008). Now, 
as all of these components contribute to the gravitational field of a 
galaxy, obtaining proper models that include the effects of all parts 
is a rather difficult problem. However, the contribution of all parts 
is limited to certain distance scales, so that not all components have 
to be included in a reasonably realistic model. Therefore, it is com- 
monly accepted that many of the main aspects of the galactic dy- 
namics can be described, in a quite approximate way, with models 
that only consider the thin galactic disc. 

Accordingly, the study of the gravitational potential gener- 
ated by an idealized thin disc is a problem of great astrophys- 
ical relevance and so, through the years, different approaches 
have been used to obtain such kind of thin disc models (see 
iBinnev & Tremainel d2008h and references therein). So, once an 
expression for the gravitational potential has been derived, corre- 
sponding expressions for the surface mass density of the disc and 
for the circular velocity of the disc particles can be obtained. Then, 
if the expression for the circular velocity can be adjusted to fit the 
observational data of the rotation curve of a particular galaxy, the 
total mass can be obtained by integrating the corresponding surface 
mass density. 
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However, although most of these thin disc models have sur- 
face densities and rotation curves with remarkable properties, many 
of them mainly represent discs of infinite extension and thus they 
are rather poor flat galaxy models. Therefore, in order to obtain 
more realistic models of flat galaxies, is better to consider methods 
that permit the obtention of finite thin disc models. Now, a sim- 
ple method to obtain the gravitational potential, the surface density 
an d the rotation c urve of thin discs of finite radius was developed 
by iHun ter ( 1963), the simple st example of a disc obtained by this 
method being the well known iKalnai si dl972h disc. 

In a previous paper dGonzalez & Reinal I2006T) we use the 
Hunter method in order to obtain an infinite family of thin 
discs of finite radius with a well behaved surface mass den- 
sity, an infinite family of generalized Kalnajs discs. Also, 
the motion of test particles in the gravitational fields gener- 
ated by the first four members of this fam ily was studied 
dRamos-Caro. Lopez-Suspez & Gonzalez! 120081) and a new infi- 
nite family of self-consistent models was obtained as a superpo- 
sition of members belonging to th e generalized Kalnajs family 
dPedraza, Ramos-Caro & Gonzalezll2008l) . 

In [Gonzalez & Reinal d2006l) . the family of disc models 
was derived by requiring that the surface density behaves as a 
monotonously decreasing function of the radius, with a maximum 
at the center of the disc and vanishing at the edge. So, although the 
mass distribution of this family of discs present a satisfactory be- 
havior in such a way that they could be considered adequate as flat 
galaxy models, their corresponding rotation curves do not present 
a so good behavior as they do not reproduce the flat region of the 
observed rotation curve. 
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On the other hand, in Pedraza, Ra mos-Caro & Gonzalez! 

d2008l) the new family of discs was obtained by superposing the 
members of the generalized Kalnajs family in order that the result- 
ing surface density can be expressed as a well behaved function 
of the gravitational potential, in such a way that the corresponding 
distribution functions can be easily obtained. Furthermore, besides 
present a well-behaved surface density, the models also present ro- 
tation curves with a better behavior than the generalized Kalnajs 
discs and are radially stable, whereas vertically unstable. Then, 
apart of the stability problems, these discs can be considered as 
quite adequate models in order to describe satisfactorily a great va- 
riety of galaxies. 

In agreement with the above considerations, in this paper we 
explore the possibility of obtain some thin disc models in which 
the circular velocities can be adjusted very accurately to fit the ob- 
served rotation curves. In order to do this, we will consider a differ- 
ent approach in order to solve the boundary value problem defining 
the gravitational potential. So, instead of assuming a given behav- 
ior of the surface density, we will obtain from the general solution 
of the Laplace equation a quite general expression for the circu- 
lar velocity, which it is expressed as a power series expansion of 
a dimensionless radial coordinate. Then, after determining the co- 
efficients of the series expansion by means of a proper numerical 
adjust, the corresponding surface densities and all the quantities 
characterizing the kinematic behavior of the particular models can 
be obtained. Also, the total mass of the galaxies can be obtained by 
integrating their surface density. 

The paper is organized as follo ws. First, in Section[2] closely 
following iGonzalez & Rein3 d2006h . we present the general finite 
thin disc model obtained by means of the Hunter method. Then, 
in Section [3] we particularize the model by taking a finite number 
of terms in the series expansion of the general model. Also we ob- 
tain explicit expressions for the kinematical quantities characteriz- 
ing the behavior of the discs. In Section|4]the models are then fitted 
to data of the observed rotati on curves of four galaxies in the Ursa 
Major cluster, as reported in|Verheiien & Sancicil bOOll) . Finally, 
in Section[5] we discuss the obtained results. 



2 GENERAL FINITE THIN DISC MODELS 

We begin by considering the Laplace equation for an axially sym- 
metric problem, written as 
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where (R, if, z) are the usual cylindrical coordinates. We will sup- 
pose also that, besides the axial symmetry, the gravitational poten- 
tial has symmetry of reflection with respect to the plane z = 0, 



(2) 



the relation 
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Also, given z), the density E(-R) of the surface distribution 
of matter can be obtained using the Gauss law and, by using the 
equation l[3}, we obtain 
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Thus, in order to have a surface density corresponding to a finite 
disclike distribution of matter, we impose boundary conditions in 
the form 
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in such a way that the matter distribution is restricted to the disc 
z = 0,0^R^a. 

In order to properly pose the boundary value problem, we 
introduce oblate spheroidal coordinates, whose symmetry adapts 
in a natural way to the geometry of the model. These coordi- 
nates are related to the u sual cylindrical coordinates by the relation 
jMorse&Fesbachlll953l) 
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where ^ £ < 00 and — 1 ^ 77 < 1. The disc has the coordinates 
£ = 0, < ?7 2 < 1. On crossing the disc, the r\ coordinate changes 
sign but does not change in absolute value. The singular behaviour 
of this coordinate implies that an even function of r\ is a continu- 
ous function everywhere but has a discontinuous r\ derivative at the 
disc. 

In terms of the oblate spheroidal coordinates, the Laplace 
equation can be written as 
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and we need to find solutions that are even functions of 77 and with 
the boundary conditions 
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According to this, the Newtonian gravitational potential for the ex- 
terior of a finite t hin disc with an axially symmetric matter density 
can be written as (Bateman 1 944b . 



which implies that the normal derivative of the potential satisfies 
the relation 



— (R,-z) = ~ — {R,z), 
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in agreement with the attractive character of the gravitational field. 
We also assume that 9 ( J>/9z does not vanish on the plane z = 0, in 
order to have a thin distribution of matter that represents the disc. 

So, given a potential &(R, z) with the previous properties, we 
can easily obtain the circular velocity v c (R), defined as the tangen- 
tial velocity of the stars in circular orbits around the center, through 



m,v) = -£c a ««in(0ft»('7), do) 



where C*2n are arbitrary constants, P2n(j]) and <?2n(£) = 
i 2n+1 Q2n(i£,) are the usual Legendre polynomials and the Leg- 
endre functions of second kind, respectively. 

With this general solution for the gravitational potential, the 
circular velocity can be written as 



~ R 2 
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while the surface matter density is given by 

1 °° 

= ^^]TC 2n (2n + l) g2 , l+1 (0)P2„(r?), (12) 



where r\ — yl — R 2 and R — r /a. So, by integrating on the total 
area of the disc, we find the expression 

MG 



= Co, 



(13) 



which allows to compute the value of the total mass M. Accord- 
ingly, all the quantities characterizing the thin disc model are de- 
termined in terms of the set of constants C 2n , which can be deter- 
mined from the observational data corresponding to rotation curves 
for some particular galaxy. 



3 A FAMILY OF PARTICULAR MODELS 

We now explore the possibility of obtain particular thin disc models 
in which the expression for the circular velocity can be very accu- 
rately adjusted with the observed data from the rotation curve of 
a given galaxy. However, in order to do this, first the sum must 
be limited to a finite number of terms. This correspond to take 
C 2 „ = for n > m, with m a positive integer. So, after replace the 
derivatives of the Legendre polynomials, the expression i ll It can be 
cast as 



(14) 



where the A 2n constants are related with the previous constants 
C 2 „ through the relation 



C 2n = 
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for n 7^ 0, which is obtained by equaling expressions dl lb 
and \14\ a nd by usin g properties of the Legendre polynomials 
(Arfken& Weber 2005). Then, if the constants A 2n are determined 
by a fitting of the observational data of the corresponding rotation 
curve, the relation J 1 5b gives the values of the constants C 2n that 
defines the particular thin disc model through JlOb . 

As we can see, the value of Co it is not determined by ex- 
pression J 1 5b - However, it is clear from d 1 2b that the surface mass 
density diverges at the disc edge , when r\ — 0, unless that we im- 
pose the condition (Hunter 1963) 



numerically adjusted to any set of data. Accordingly, expression 
j 14b can be considered as a kind of "universal rotation curve" for 
flat galaxies, which can easily be adjusted to the observed data of 
the rotation curve of any particular spiral galaxy. 

Now then, besides the circular velocity, there are two other im- 
portant quantities concerning the kinematics of the models, which 
describe the stability against ra dial and vertical perturba tions of 
particles in quasi-circular orbits (Binnev & Tremaine 2008). These 
two quantities, which must be positive in order to have stable cir- 
cular orbits, are the epicyclic or radial frequency, defined as 
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and the vertical frequency, defined as 
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is the effective potential and £ = Rv c is the specific axial angular 
momentum. 

By using then l(4j and J20b in d 1 8b . we can easily obtain the 
relation 



k 2 (R) 



1 dvi 2vi 
R~dR + "R 2 



(21) 



in such a way that, by using J 1 4b . the epiciclic frequency can be 
cast as 
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where k = an. Also, from the expression for the Laplace operator 
in cylindrical coordinates and using 10, J19b and ( |20b , is easy to 
see that 



(23) 



so, as the potential is a solution of the Laplace equation, by using 
j 14b the vertical frequency can be written as 
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where v = av. 
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that, after use the properties of the Legendre functions, leads to the 
expression 



c o = E(~ 1 )" +1C2 " 



(17) 



which gives the value of Co, and then of the total mass M, in terms 
of the A 2n . 

The previous expressions imply then that any particular thin 
disc model will be completely determined by a set of constants 
A 2 „, which must be choosed in such a way that the circular ve- 
locities can be adjusted very accurately to fit the observed rotation 
curves. Furthermore, as the powers R 2n are a set of linearly in- 
dependents functions, the expression j 14b is quite adequate to be 



4 FITTING OF DATA TO THE MODELS 

In order to adjust the previous model to real observed data, we 
choose four spiral galaxies of the Ursa Major cluster, the galaxies 
NGC3877, NGC3917, NGC3949 and NGC4 010. The correspond- 
ing da ta are taken from the recent paper by IVerheiien & Sancicil 
(2001), which presents the results of an extensive 21 cm-line syn- 
thesis imaging survey of 43 galaxies in the nearby Ursa Major clus- 
ter using the Westerbork Synthesis Radio Telescope. For each rota- 
tion curve data, we take the value of a as given by the last tabulated 
value of the radius. Accordingly, we are assuming that the radii of 
the galaxies are defined by the last observed data. Then, by taking 
the radii normalized in units of a, we make a non-linear least square 
fit of the data with the general relation J 1 4b . by considering in each 
case a value of m less than the number of available data points. 
In FigureQ]we show the adjusted rotation curves for the four 
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Table 1. Constants C 2 n [km 2 s 2 ] and values of m. 



NGC3877 





NGC3877 


NGC3917 


NGC3949 


NGC4010 


m 


6 


7 


5 


7 


C 


17452.78 


11258.92 


15686.77 


8801.12 


C 2 


26564.93 


17210.85 


24051.95 


12875.61 


c 4 


13926.79 


8998.92 


13011.66 


4944.66 


c 6 


7478.82 


4573.34 


7127.38 


1117.08 


C 8 


4053.13 


2213.64 


4577.68 


866.64 


Cio 


1887.23 


1063.01 


2096.77 


1019.52 


C12 


498.29 


640.40 




818.98 


Cl4 




264.69 




419.19 



Table 2. Morphological type, radius a and total mass M . 





Type 


a [kpc] 


M [kg] 


M [M e ] 


NGC3877 


Sc 


11.74 


9.47 x 10 40 


4.76 x 10 10 


NGC3917 


Scd 


15.28 


7.95 x 10 40 


3.95 x 10 10 


NGC3949 


Sbc 


8.72 


6.32 x 10 40 


3.18 x 10 10 


NGC4010 


SBd 


10.84 


4.41 x 10 40 


2.22 x 10 10 



galaxies considered. The points with error bars are the observa- 
tions, as reported in IVerheiien & Sancicil fcOOlh . while the solid 
line is the rotation curve determined from d 1 4b with the values for 
the A 2n given by the numerical fit. As we can see, the relation d 14b 
fit quite accurately to the observed data of the four galaxies consid- 
ered. Then, from the obtained values for A 2n , the corresponding 
values of the C 2n are determined by using d 1 5b and d 1 7b - In Ta- 
ble [T] we present the values of C 2n for the four galaxies as well as 
the corresponding value of m used in d 14b - In Tablef2]we indicate, 
for each galaxy, the morphological type according to the Hubble's 
classification of galaxies, the radius a in kpc and the total mass M, 
both in kg and in solar mass units (.Mo). 

Now, as the set of constants C 2n it defines completely each 
particular thin disc model, we can easily compute all the physical 
quantities characterizing each galaxy. However, as explicit expres- 
sions for the gravitational potential &(R, z) and the surface mass 
density E(i?) can be easily obtained by using the values of the C 2n 
at expressions d 1 0b and d 1 2b . we will not present them here. Insteed, 
we plot in Figure [2] the surface densities for the four galaxies, as 
functions of the dimensionless radial coordinate R = R/a. For the 
four galaxies we obtain a well behaved surface mass density, hav- 
ing a maximum at the disc center and then decreasing until vanish 
at the disc edge. 

In a similar way, we can compute the epiciclic and vertical 
frequencies by using ( [22b . l !24b and the values of the constants A 2n 
obtained from the numerical fit. However, as with the surface mass 
densities, we will not present the explicit expressions here and, in- 
stead, we only show the corresponding plots. So, in Figure [3] we 
show the plots of the epiciclic frequencies for the four galaxies 
considered and, in Figure [4] the corresponding plots of the verti- 
cal frequencies. From the plots at Figure[3]we can see that only the 
galaxy NGC4010 presents a small region of radial instability near 
the disc edge. On the other hand, as it is shown at Figure[4] the four 
galaxies are instable against vertical perturbations. 




0.4 0.6 

NGC3949 




Figure 1. Plots of the circular velocities v c in km/s, as functions of the di- 
mensionless radial coordinate R = R/a, for the spiral galaxies NGC3877, 
NGC3917, NGC3949 and NGC4010. 
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NGC3877 NGC3877 




Figure 2. Plots of the surface densities £ X 10 3 in kg/m 2 , as functions Figure 3. Plots of the epiciclic frequencies k 2 x 10 5 in (km/s) 2 , as 

of the dimensionless radial coordinate R = R/a, for the spiral galaxies functions of the dimensionless radial coordinate R = R/a, for the spiral 

NGC3877, NGC3917, NGC3949 and NGC4010. galaxies NGC3877, NGC3917, NGC3949 and NGC4010. 
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NGC3877 




0.2 0.4 0.6 0.8 1 



NGC3917 




0.2 0.4 0.6 0.8 1 



NGC3949 




0.2 0.4 0.6 0.8 1 



NGC4010 




0.2 0.4 0.6 0.8 1 



Figure 4. Plots of the vertical frequencies v 2 X 10 4 in (km/s) 2 , as func- 
tions of the dimensionless radial coordinate R = R/a, for the spiral galax- 
ies NGC3877, NGC3917, NGC3949 and NGC4010. 



5 DISCUSSION 

We presented four particular thin disc models adjusted in order to 
accurately fit the observed data of the rotation curves for the galax- 
ies NGC3877, NGC3917, NGC3949 and NGC4010 of the Ursa 
Major cluster. These models present well behaved surface densi- 
ties, that resembles the observed luminosity profile of many spiral 
galaxies, and the obtained values for the corresponding total mass 
Ai it agrees with the expected order of magnitude. Accordingly, the 
here obtained expression for the circular velocity, equation d 14b . 
can be considered as a kind of "universal rotation curve" for flat 
galaxies, which can easily be adjusted to the observed data of the 
rotation curve of any particular spiral galaxy. 

On the other hand, in one of the models we obtain a small 
region near the disc edge with instability agains radial perturba- 
tions. Now, as the models are completely determined by the set of 
constants A2 n , which are fixed by the numerical fit of the rotation 
curve data, there are not free parameters that can be adjusted by 
requiring radial stability. A posible solution for this problem may 
be consider a less restrictive numerical fit that leaves some free pa- 
rameters. This can be done by taking the summation in expression 
J 14b until a value of m greater than the available number of data 
points. 

However, the models present a central region with strong 
instabilities against vertical perturbations of particles in quasi- 
circular orbits. This result was expected as a consequence of the 
fact that the models only consider the thin galactic disc. Indeed, as 
we can see from expression j23\ , vertical instability will be always 
present in models constructed from solutions of Laplace equation 
and adjusted in such a way that their circular velocities reproduce 
the central region of the observed rotation curves, where the veloc- 
ity rises linearly with the radius. Therefore, more realistic models 
must be considered that include the non-thin character of the galac- 
tic disc, or the mass contribution of the spheroidal halo. 

In agreement with the above considerations, we can consider 
the set of models here presented as a first approximation to the 
obtaining of quite realistic models of spiral galaxies. In particular, 
we believe that the values of A4 that were obtained for the four 
galaxies studied may be taken as a quite accurately estimative of 
the mass upper bound of these galaxies, since in the model was 
considered that all their mass was concentrated at the galactic disc. 
Accordingly, we are working now in a generalization of the model 
that includes the mass contribution of the spheroidal halo, in such 
a way that we can overcome the vertical instability problem and 
obtain some estimative of the relative contributions of the halo and 
the disc to the total mass of the galaxies. 
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